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VELOCITY  SHEAR  STABILIZATION  OF  THE  CURRENT  CONVECTIVE  INSTABILITY 

I.  INTRODUCTION 

Large  scale  ftructure  in  plasma  density  enhancements  In  the  auroral  F 
layer  Ionosphere  has  been  a  topic  of  continued  research  over  the  past 
several  years  (Vickrey  et  al.,  1980;  Keskinen  and  Ossakow,  1982;  1983a, b; 

Tsunoda  and  Vickrey,  1983).  Theoretically  these  so  called  F  region  plasma 
"blobs"  have  been  shown  to  be  susceptible  to  the  E_  x  B_  gradient  drift 

instability  (EBI)  (Linson  and  Workman,  1970)  and  in  an  _E_  x  stable 

geometry  to  the  current  convective  instability  (CCI)  (Ossakow  and 
Chaturvedi,  1979),  or  a  combination  of  the  two  (Vickrey  et  al. ,  1980; 

Keskinen  and  Ossakow,  1982).  The  current  convective  instability  was 
originally  dlcussed  in  the  context  of  the  positive  column  of  a  laboratory 
gas  discharge  by  Lehnert  (1958),  and  Kadomtsev  and  Nedospasov,  (1960). 
Ossakow  and  Chaturvedi  (1979)  suggested  that  the  diffuse  auroral 

A 

precipitation,  constituting  a  weak  magnetic  field  aligned  current  (z), 
could  destabilize  an  inhomogeneous  plasma  having  the  density  gradient 

A 

pointing  in  the  direction  perpendicular  (northward,  x)  to  the  magnetic 
field  .  The  unstable  waves  have  wavevectors  pointing  in  the  direction 
perpendicular  (eastward,  y)  to  both  the  ambient  geomagnetic  field  and  the 
density  gradient.  DNA  Wideband  satellite  observations  (Rlno  et  al. ,  1978) 
indicated  that  a  regularly  occurring  scintillation  enhancement  due  to  F 
region  ionospheric  irregularities  can  be  identified  in  the  nighttime 
auroral  zone  data.  This  data  is  characterized  by  diffuse  auroral  particle 
precipitation  and  total  electron  content  (TEC)  gradient  pointing 
northward.  Furthermore,  the  dominant  modes  appear  to  be  in  the  north-south 
direction  parallel  to  the  density  gradient  in  contrast  to  the  predictions 
of  the  linear  theory.  Indeed,  Chaturvedi  and  Ossakow  (1979)  have  shown 
that  waves  parallel  to  the  density  gradient  are  generated  by  a  nonlinear 
mode  coupling  mechanism. 

HuMMript  approved  November  1,  1983. 
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There  is  increasing  evidence  chat  a  strong  velocity  shear  exists  in 
the  diffuse  auroral  region  (Kelley  et  al.,  1976;  Fejer  et  al.,  1983;  S. 
Basu,  private  communication,  1983).  The  east-west  plasma  velocity  is 
sheared  as  one  moves  in  the  north-south  direction.  The  coupling  to  the 
neutral  atmosphere,  where  the  neutral  flow  velocity  usually  is  strongly 
sheared,  probably  generates  this  plasma  velocity  shear.  Recent 
calculations  have  shown  that  velocity  shear  has  .  a  strong  stabilizing 
influence  on  interchange  instabilities,  such  as  the  Rayleigh-Taylor 
instability  (Hamieri,  1979;  Guzdar  et  al.,  1982;  1983)  and  the  x  B_ 
instability  (Perkins  and  Doles,  1975;  Huba  et  al. ,  1983)  leading  one  to 
anticipate  similar  effects  in  the  case  of  the  CCI  also. 

In  this  paper  we  investigate  the  effects  of  velocity  shear  on  the  CCI 
and  its  role  in  producing  large  scale  east-west  structure  in  ionospheric  F 
layer  plasma  blobs  (Tsunoda  and  Vickrey,  1983).  We  find  that  velocity 
shear  stabilizes  the  short  wavelength  modes  (k^L  >  1,  where  ky  is  the  east- 
west  mode  number,  L  is  the  density  gradient  scale  length)  and  thus  excites 
a  long  wavelength  mode. 

The  organization  of  this  paper  is  as  follows.  In  the  next  section  we 
present  the  general  nonlocal  theory  and  present  its  results  in  the  third 
section.  We  discuss  the  results  and  apply  them  to  the  diffuse  auroral 
ionosphere  in  the  fourth  section,  and  the  fifth  section  contains  the 
conclusions. 


II.  THEORY 


The  geometry  we  use  is  as  follows.  The  background  magnetic  field  and 

A  A 

the  field  aligned  currents  are  along  the  z-direction  (B.  *  B  £;  *  V^jb). 

The  density  gradient  and  the  inhomogeneous  electric  field  [Eox(*)l  are 
along  the  x-direction  (northward).  The  Inhomogeneous  electric  field  gives 
rise  to  an  E_ x  ^  sheared  flow  along  the  negative  y  direction  (eastward). 

The  basic  assumptions  used  in  the  analysis  are  as  follows:  (i)  the 

perturbed  quantities  vary  as  fif  -  <5f(x)  exp[i(k^y  +  k^z  -  a»t)],  where  ky 

and  k2  are  the  wavenumbers  along  the  y  nd  z  directions  and  to  “  a>r  +  iy  , 

implying  growth  for  y  >  0  ;  (ii)  the  ordering  in  the  frequencies  is  such 

that  to,  v£n(vei)  <<  Q.  (8@)«  where  v^n  and  v  A  are  the  ion-neutral  and 

electron-ion  collision  frequencies  and  ft  and  ft,  are  the  electron  and  ion 

e  i 

cyclotron  frequencies,  respectively;  (iii)  we  retain  the  ion  inertia 
effects,  thereby  including  the  ion  polarization  drift,  but  ignore  the 
electron  inettla;  and  (iv)  we  ignore  f inite-Larmor  radius  effects  by 
limiting  the  wavelength  domain  to  kpt  «  1,  where  is  the  mean  ion  Larmor 
radius. 

A  key  feature  of  our  analysis  is  that  a  nonlocal  theory  is 
developed.  That  is,  the  mode  structure  cc  the  potential  in  the  x 
direction,  the  direction  in  which  the  density  and  flow  velocity  are  assumed 
to  vary,  is  determined  by  a  differential  equation  rather  than  an  algebraic 
equation  obtained  by  Fourier  analysis.  This  technique  allows  one  to  study 
modes  which  have  wavelengths  comparable  to  the  scale  sizes  of  the 
inhomogeneities  (i.e.,  kyL  <  1,  where  L  represents  the  scale  lengths  of  the 
boundary  layer).  In  fact,  nonlocal  theory  is  essential  to  describe  the 
Kelvin-Helmholtz  instability  driven  by  a  transverse  velocity  shear 
(Mikhailovskii,  1974). 

Based  on  the  assumptions  described  above,  the  fundamental  fluid 
equations  used  in  the  analysis  are  continuity  and  momentum  transfer  in  the 
neutral  frame  of  reference 


0  -  -e(E  +  (l/c)V^  x  B )  -  «v#iV€ 
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(2) 
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where  m(M)  is  the  electron  (ion)  mass,  subscript  j  denotes  species  (e  for 
electrons  and  i  for  ions),  V  is  velocity  and  E_  is  electric  field,  n  is 
density,  c  is  speed  of  light,  and  e  is  the  absolute  value  of  the  electron 
charge.  From  eqs.  (1)  and  (2)  we  obtain  the  zeroth  order  perpendicular 
electron  and  ion  velocities  as 


V  -  +  (c/B)E  x  z  +  0(v  ./n  )  ,  (4) 

—el  —  —  ei  e 

-ii“  +<c/B  [li  *  1  +  Kn/fli)ii  +  (1/Qi)  (VV)  Ei'  *  (5) 

where  eB/Mc  is  the  ion  cyclotron  frequency. 
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Taking  the  dot  product  of  eqs.  (1)  and  (2)  with  z_  ,  we  obtain  the 
zeroth  order  parallel  velocities  as 


•*  Eoz/avei’ 


'  «Eoz/Mvin‘ 


(6) 

(7) 


We  perturb  equations  (l)-(3)  and  linearize  them  in  the  electrostatic 

limit.  We  substitute  E  •  -V<4  ■  -  V($  (x)  +  5$),  V  -V  +  6V  ,  and 

—  o  — a  — r*  —a 

n  •  nQ(x)  +  6n  into  eqs.  (l)-(3).  Since  the  equilibrium  quantities  vary 

in  the  x  direction  we  do  not  Fourier  analyze  in  the  x  direction.  We 

Fourier  analyze  in  the  y  and  z  directions.  Using  eqs.  (4)  and  (5)  along 
with  eqs.  (1)  and  (2)  we  obtain  the  perturbed  velocities  in  the 

perpendicular  direction  as 
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5Xsi"  +(c/B)  5ex  x  z  +  o(vel/ne>  , 
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sviy  -  (C/3)  [£■  0> - 2L2Z - ±2_  kv5^  ] 


•3_ 

3x 


Q, 
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where  VQy  ■  -(c/B)  EQX.  From  eqs.  (1,2)  and  (6,7)  the  perturbed  parallel 
velocities  are  given  as 


SV 


ez 
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-e  5Ez/mvel, 
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z  m 


(10) 

(11) 


We  substitute  eqs.  (4)  thru  (11)  in  eq.  (3)  and  assume  quasi- 

neutrality,  <5n  «  5n  «  <5n  ,  to  obtain  the  mode  structure  equation  for  the 

G  X 

perturbed  electrostatic  potential  as  (where  we  have  dropped  6  from  5$  and 
indicate  the  perturbed  potential  by  $  ) 


+  p$'  +  q$  ■  0 


where 


P  -  t  +  U»la/»j)  k2v0,/iil 


(12) 


(13) 
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where  we  assumed  V"  /a.  «  1  and  ignored  terms  of  0(v  . /n  )  and 
oy  i  ei  e 


V  -  (V  -V.  ),  e  -  n'/n  , 
oz  ez  iz  o  o  * 


I*  *  ik|  (a./o.i)  -  y 


(15) 


y.  ■  «  -  k  V  -  1c  V  ,  and  au  „  ■  ~  v.  . 

y  oy  z  oz  *  2  in 


a>  +  iv. 


The  primes  indicate  derivatives  with  respect  to  x.  Equation  (12)  is 
similar  to  the  equations  that  have  been  derived  in  the  context  of  the 
Rayleigh-Taylor  (Guzdar  et  al. ,  1982;  1983)  and  E_ x  instabilities 
(Huba  et  al. ,  1983).  In  the  absence  of  the  velocity  shear 
(V'y  ■  V''  ■  0)  the  first  three  terms  in  q(x)  (eq.  14)  yield  the  nonlocal 
current  convective  instability.  Equation  (12)  is  solved  numerically  in 
various  parameter  domains  and  the  results  are  presented  in  the  next 
section. 
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III.  ANALYSIS 


Before  solving  eq.  (12)  we  consider  some  known  limiting  cases.  First, 
in  the  absence  of  field  aligned  currents  (VQ2  *  0),  for  a  homogeneous 
(7n  ■  0)  collisionless  (vin  ■  0)  plasma  eq.  (12)  yields  the  well  known 
transverse  (k2  *  0)  Kelvin-Helmholtz  instability  (Mikhailovskii,  1974) 

♦"  +  (-  +  (.  -Vv  )  )»  ■  0  (16) 

y  °y 

The  Kelvin-Helmholtz  instability  is  shown  in  Fig.  (1)  curve  A  (Michalke, 
1964).  The  velocity  profile  used  is  VQy  ■  VQy  tanh(x/L).  We  see  from  the 
figure  that  the  collisionless  Kelvin-Helmholtz  mode  is  unstable  in  the 
domain  0  <  kyL  <  1  with  the  growth  rate  maximizing  at  kyL  -  1.5  with  a 

maximum  value  of  0.18  V  /L. 

oy 

In  the  absence  of  the  velocity  shear  we  obtain 


+  £i 


+  r-  k2  +  ik2 

L  y  z  a. 


k  v  n,  n 

+  — - — -  (ik2  —  -  k 

u>  ui+iv,  ^  z  v  y  ' 
in  e 


(17) 


Equation  (17)  describes  the  nonlocal  current  convective  instability.  In 
the  local  limit,  3/3x  «  ky,  we  obtain  the  local  dispersion  equation 
(Chaturvedi  and  Ossakow,  1981) 


u.2  +  u[ivin  +  i(k2/lt2)(Qe01/vei)]  -  (k2/ky)  -  0.  (18) 

The  solution  of  eq.  (18)  yields  the  growth  rate  (Chaturvedi  and  Ossakow, 
1981)  with  u  *  a>r  +  iy  as 


y 


>-kz/ky)  ^ei/ne>  Voze 
(kz2/ky)  +  (veini/vin^ 


(19) 


for  vt*  >  4  (kx/ky)  Voxe  where  viQ-  vit  +  (k2/kj)  •  Note 
that  for  instability  k2/ky,  vei,  V  qz,  and  e  have  to  be  nonzero.  In 
fact,  the  optimum  value  of  k2/ky  can  be  calculated  by  maximizing  eq.  (19) 
with  respect  to  k2/ky. 
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The  growth  rate  maximizes  at 


kz  ■  k,  (vvtn/a.ai)1/2- 

with  a  maximum  value 

V*  •TVoze(',eiai/vln  ae>  ^ 

For  v  ./Q  »  10-4,  v.  /n.  -  10  4, 
ex  e  in  x 


(20) 


(21) 


(y/V  e)  -  0.5 
k 1  oz  Jmax 


(22) 


Now  we  solve  eq.  (17)  numerically  for  various  k  "  ising  the  density 


profile 


nQ(x)  -  nQ(l  +  6  tanh(x/L) )/ (1-6)  (23) 

shown  In  figure  2  for  6  -0.8. 

Figure  (1)  curve  B  shows  the  growth  rate  normalized  to  V  /L  as  a 
function  of  kyL  for  6-0.8  and  kz  given  by  eq.  (20).  We  note  that  the 
growth  rate  increases  monotonically  with  kyL  and  asymptotes  to  the  local 
growth  rate  (eq.  22). 

We  now  introduce  the  velocity  shear  and  solve  eq.  (12)  numerically. 
We  use  the  density  profile  given  in  eq.  (23).  The  velocity  profile  is 
determined  self-cons latently  from  the  zeroth  order  continuity  equation 

-  0  •  (24) 

The  eigenvalues  are  determined  for  various  values  of  the  parameter 
* 

s(s  Voy/Vo8),  which  serves  as  a  measure  of  the  velocity  shear.  In  Fig.  (3) 

we  plot  the  growth  rates  normalized  to  the  shear-free  current  convective 

growth  rate,  y/(V  /L),  as  a  function  of  the  normalized  wavenumbers, 
08  -4  -4 

kyL.  We  choose  *10  ,  v^/Q^  -  10  ,  and  6  -  0.8.  Curve  A  gives 

the  case  s  •  0  (same  as  curve  B,  Fig.  (1)).  Curves  B,  C  and  D  are  for  s  - 
0.S,  1,  and  2,  respectively.  We  can  draw  several  conclusions  from  this 
figure.  First,  velocity  sheer  stabilizes  the  short  wavelength  modes 

8 


(Perkins  and  Doles,  1975;  Guzdar  et  al.  ,  1982;  1983;  Huba  et  al.  ,  1983). 

Second,  the  growth  rate  maximizes  at  k^L  ~  1.5  for  s  *  0.5;  and  as  s  is 

increased  the  peak  growth  rate  decreases  and  the  growth  rate  maximizes  at 

longer  wavelengths.  For  s  *  2,  the  growth  rate  maximizes  at  kyL  -  0.55 

with  y  ~  0.083(V  /L).  Third,  when  s  ■  1,  i.e.,  when  the  horizontal 

flow  velocity  is  of  the  order  of  the  vertical  plasma  drift  velocity,  modes 

with  X  ~  (2it/0.9)L  are  the  fastest  growing  modes  (v  “  0.142 

'max 

Vrt_/L  -  0.142  Vrt„/L) .  We  note  that  as  s  +  ®  the  V'  and  V"  terms 
oz  oy  o  o 

dominate  in  eq.  (14)  and  the  growth  rate  curve  peaks  around  k^L  ~  0.45  (not 
shown  in  the  figure)  resembling  curve  A  of  fig.  1. 


•TfW- 


IV.  DISCUSSION 


Recent  experiments  conducted  at  Chatanika  (Tsunoda  and  Vickrey,  1983) 
indicate  that  high  latitude  F-region  plasma  blobs  have  an  east-west  plasma 
density  structure  with  scale  sizes  of  130  km  and  greater.  Even  though 
these  experiments  did  not  measure  the  velocity  shear  in  this  region,  other 
evidence  (Kelley  and  Carlson,  1977)  indicates  the  presence  of  a  strong 
east-west  flow  velocity  shear  of  the  order  1-20  sec”*  in  the  diffuse 
auroral  region.  Assuming  that  the  east-west  flow  is  ~  250  m/s  and  the 
shear  scale  length  (of  the  order  of  the  density  gradient  scale  length  jLjj  * 
Lv)  is  25  km.  Fig.  (3)  shows  that  the  generalized  current  convective 

A 

instability  has  perpendicular  scale  sizes  of  order  150  km  (s  ■  1;  VQy  * 
Voz;  kyL  *  0.9).  The  corresponding  growth  rate  is  given  as 
Y  -  0.142  (V^/L)  "*  0.142  10  ^s  *  (around  10  minutes  e-folding  time).  For 
larger  shears,  for  example,  with  the  east-west  flow  velocity,  VQy,  twice 
that  of  the  field  aligned  flow,  VQ2,  the  instability  generates 
irregularities  of  scale  sizes  around  300  km,  but  with  a  smaller  growth 
rate,  0.08  (VQZ/L)  (20  min  e-folding  time  scale  for  VQZ/L  around  10“  s-*). 

Since  velocity  shear  has  similar  effects  on  the  E_  x  B_  gradient  drift 
instability  (Perkins  and  Doles,  1975;  Huba  et  al. ,  1983)  it  is  important  to 
compare  our  present  results  with  the  results  of  Huba  et  al.  (1983). 
Referring  to  Fig.  (5)  of  the  paper  by  Huba  et  al. ,  we  see  that  east-west 
flow  comparable  to  the  equilibrium  north-south  flow  (due  to  the  horizontal 
electric  fields  generating  the  _E_  x  B_  instability)  excites  modes  with 
wavelengths  -  60  km  (kyL  *  2.6,  L  ■  25  km).  The  growth  rate  of  these  modes 
is  ~  0.5  (Vox/L)  with  e-folding  times  of  the  order  3  minutes  (for  Vox/L 
10  *).  From  this  we  can  see  that  in  the  diffuse  auroral  region,  the 
generalized  CCI  and  EBI  (CCI  and  EBI  with  velocity  shear)  can  excite  long 
wavelength  modes  with  east-west  scale  sizes  of  100  km  or  greater.  However, 
we  point  out  that  in  the  absence  of  the  E_  x  J3_  instability  (zero  westward 
electric  field  in  the  presense  of  say  a  poleward  density  gradient  on  a 
blob)  the  soft  particle  precipitation  causing  in  turn  currents  of  moderate 
strength  ■-  250  m/s  (4  yA/m2,  n  ■  10**  m”^)  can  still  excite  the  long 
wavelength  modes. 
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V.  CONCLUSIONS 


We  have  studied  the  effects  of  velocity  shear  on  the  current 
convective  instability.  We  have  shown  that  velocity  shear  suppresses  the 
short  wavelength  modes  and  preferentially  excites  a  long  wavelength  mode 
with  kyL  <  1.  This  suggests  that  a  generalized  current  convective 

instability  (CCI  with  east-west  flow  velocity  shear)  with  currents 
~  4yA/m2  ('•  250  m/s  parallel  flow)  can  possibly  excite  modes  with 

wavelengths  of  150  km  or  more.  Since  the  x  _B_  instability  can  also  excite 
long  wavelength  modes(A  -  60  km  with  a  westward  electric  field  ~  lOmV/m, 
and  north-south  flow  ~  east-west  flow),  we  conjecture  that  the  soft 
particle  precipitation  alone  can  destabilize  a  stable  _E_  x  ^  diffuse  auroral 
F  region  geometry  by  the  CCI  (Ossakow  and  Chaturvedi,  1979).  In  turn.  The 
CCI  in  the  presence  of  velocity  shear  can  generate  irregularities  in  the 
long  wavelength  domain. 

We  note  that  a  unified  treatment  of  current  convective  and  E_  x  _B 
instabilities,  with  velocity  shear  and  magnetic  shear  and  nonlinear  effects 
included,  is  needed  in  order  to  come  to  a  better  understanding  of  the 
morphology  of  the  irregularities  evident  in  auroral  F  region  plasma 
enhancements.  This  subject  is  presently  under  investigation  and  will  be 
presented  in  a  future  paper. 
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Figure  1 


.  Normalized  growth  rate  versus  kyL  curves  for  Relvin-Helmholtz 
(curve  A)  and  current  convective  (curve  B)  instabilities.  Curve 


A  refers  to  velocity  profile  V- 


tanh  x/L.  Curve  B 


oy  oy 

refers  to  the  density  profile  of  eq.  (23),  with  6  ■  0.8,  and 


Vin/S1i 


10 


-4 


vei/ne 


10 


-4 
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Generalized  current  convective  instability  as  a  function  of 

A  a 

s(=  V  /V  ).  Curve  A  refers  to  s  ■  0.  Curves  B,  C,  and  D 
v  oy  oz ; 

refer  to  s  ■  0.5,  1,  and  2,  respectively.  We  used  6  *  0.8, 
v./Q  *  10  \  and  v,  /q,  ■  10 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OICY  ATTN  I.  GARFUNKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINCTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  JUTLAND  DRIVE 
AUSTIN,  TX  78758 
OICY  ATTN  L.  SLOAN 
OICY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON.  TX  75080 
OICY  ATTN  H.  LOCSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  WUshira  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  ARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 
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general  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
Fairbanks,  ak  99701 
(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
OICY  ATTN  T.M.  DAVIS  (UNCLASS  ONLY) 
OICY  ATTN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-SASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAICN,  II  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  2231 1 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTI  TEL  &  TELEGRAPH  CORPORATION 

500  Washington  avenue 

NUTLEY,  SJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

UOa  TORRIYANA  ROAD 
P.O.  BOX  85154 
SAN  DtECO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIES  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

Q1CY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATS  STREET  (P.O  DRAWER  QQ/ 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS  I  AC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O,  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVtD  M.  TOWLE 
OICY  ATTN  l.  LOUCHLIN 
OICY  ATTN  0.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
Old  ATTS  J.  HOULE 
Old  ATTN  GEORGE  NROZ 
OICY  ATTN  W.  OLSON 
Old  ATTN  R.U.  HALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
Old  ATTN  W.F.  CREVIER 
01CV  ATTN  STEVEN  L.  CUTSCHE 
Oicr  ATTN  R.  BOGUS  CU 
Old  ATTN  R.  HENDRICK 
Old  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
Old  ATTN  F.  FACES 
Old  ATTN  H.  SCHEIBE 
Old  ATTN  CONRAD  L.  LOMGMIRE 
Old  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  87106 
Old  R.  STELLINGWERF 
Old  M.  ALME 
Old  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

Old  ATTN  JOHN  MORGANSTERN 
Old  ATTN  G.  HARDING 
OiCY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN..  VA  22101 
Old  ATTN  W.  HALL 
Old  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTS  IONOSPHERIC  RESEARCH  LAE 

PHOTOMETRICS ,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  K0FSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 
OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
Old  ATTN  FORREST  CILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  RARE  AS 
Old  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 
OICY  ATTN  L.  DeRAND 
OICY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

52S  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS ,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICT  ATTN  E.  FRIEMAN 
OICT  ATTN  E.A.  STRAKER 
OICT  ATTN  CURTIS  A.  SMITH 
OICT  ATTN  JACK  NCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  9402S 

OICT  ATTN  DONALD  NEILSON 
OICT  ATTN  ALAN  BURNS 
OICT  ATTN  G.  SMITH 
OICT  ATTN  R.  TSUNODA 
OICT  ATTN  DAVID  A.  JOHNSON 
OICT  ATTN  WALTER  G.  CHESNUT 
OICT  ATTN  CHARLES  L,  RINO 
OICT  ATTN  WALTER  JAYE 
OICT  ATTN  J.  VICKREY 
OICT  ATTN  RAT  L.  LEADABRAND 
OICT  ATTN  G.  CARPENTER 
OICT  ATTN  G.  PRICE 
OICT  ATTN  R.  LIVINGSTON 
OICT  ATTN  V.  GONZALES 
01CY  ATT.V  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WICCINS  AVENUE 
BEDFORD.  MA  01730 

OICT  ATTN  W.P.  SOqUIST 

TOTON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA.  CA  93111 
OICT  ATTN  JOHN  ISE,  JR. 

OICT  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  STS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICT  ATTN  R.  K.  PLEBOCH 
OICT  ATTN  S.  ALTSCHULER 
OICT  ATTN  D.  DEE 
OICT  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICT  ATTN  W,  REIDT 
OICT  ATTN  J.  CARPENTER 
OICT  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELINC  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OP  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  2037S 

Dr.  P.  MANGE  -  CODE  4101 
Dr.  P.  GOODMAN  -  CODE  4180 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  U.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  0221S 
DR.  J.  AARONS 

C0R.7ELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAI3RIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY.  NWl 6-262 
CAMBRIDGE,  MA  02139 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUBIN 

DR.  N.  MAYNARD  -  CODE  696 
COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  921S2 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 
ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  C.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S-  NISBET 
DR.  P.R.  R0HR3AUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVAKS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  E.  FRIEMAN 
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STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

OR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
DR.  L.E.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NN  37545 
DR.  M.  FONGRATZ 
DR.  D.  SIMONS 
DR.  G.  BARASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 

UNIVERSITY  OF  MARYLAND 
COLLECE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  GRESNWALD 
DR.  C.  MESG 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 


UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  RESEARCH  SCIENCES 
P.O.  BOX  688 
RICHARDSON,  TX  75080 
DR.  R.  HEEL IS 
DR.  V.  HANSON 
DR.  J.P.  McCLURE 

UTAH  STATE  UNIVERSITY 
4TH  AND  8T8  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  E.  BARER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
PLASMA  PHYSICS  PROGRAMME 
AHMSDABAD  380  009 
INDIA 

P.J.  PATHAX,  LIBRARIAN 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGY  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
J3AN  VARTAN  HELLMAN, 
REFERENCE  LIBRARIAN 
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